Attorney Docket Number NC17122 

Patent Application Papers of KADOUS, Tamer 



Express Mail No.: EL664600406US 



CHANNEL ESTIMATION IN A COMMUNICATION SYSTEM 

This application claims the benefit of U.S. Provisional 
Application No. 60/171,470, filed December 22, 1999. 

FIELD OF THE INVENTION 

The present invention relates generally to methods and 
apparatus for estimating a channel susceptible to distortion in a 
communication system. More particularly, the present invention relates to 
an apparatus and an associated method, for estimating channels in 
orthogonal frequency division multiplexed (OFDM) communication 
systems. 

BACKGROUND OF THE INVENTION 

Digital communication techniques have been developed and 
implemented in communication systems, including communication systems 
utilizing radio channels. Digital communication techniques generally 
permit the communication system in which the techniques are 
implemented to achiever greater transmission capacity as contrasted to 
the capacity available with conventional analog communication 
techniques. 

A communication system generally comprises a sending station 
and a receiving station communicating by way of one or more 
communication channels. Data to be communicated by the sending 
station to the receiving station is converted, if necessary, into a form to 
permit its transmission on the communication channel. A communication 
system can be defined by almost any combination of sending and 
receiving stations, including, for instance, circuit board-positioned sending 
and receiving elements as well as more conventionally-defined 
communication systems including users spaced at great distances apart 
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communicating data between each other by transmission over radio 
channels. 

When data transmitted on a communication channel is received 
at the receiving station, the receiving station acts upon, if necessary, the 
5 received data to recreate the informational content of the transmitted data. 
In an ideal communication system the data received at the receiving 
station is identical to the data transmitted by the sending station. 
However, in reality, much of the data may be distorted during its 
transmission on the communication channel. Such distortion distorts the 
10 data as received at the receiving station. If the distortion is significant, the 
informational content of portions of the data may not be recoverable. 

Q A radio communication system is one example of a 

7^ communication system utilized to transmit data between sending and 
t=J receiving stations. In a radio communication system, the communication 
m channel is formed of a radio communication channel. A radio 
communication channel may be defined within a portion of the 
fU electromagnetic spectrum. In a wireline communication system, in 
il- contrast, a physical connection between the sending and receiving 
iJf stations is implemented to form the communication channel. Transmission 
20 of data upon a radio communication channel is particularly susceptible to 
distortion, due in part to the propagation characteristics of the radio 
communication channel. Data communicated on conventional wireline 
channels are also, however, susceptible to distortion in manners 
analogous to the manner by which distortion is introduced upon the data 
25 transmitted in a radio communication system. 

In a communication system, which utilizes digital 
communication techniques, information, which is to be communicated, is 
digitized to form digital bits. The digital bits are typically formatted 
according to a formatting scheme. Groups of the digital bits, for example, 
30 are assembled to form a packet of data. 
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Orthogonal Frequency Division Multiplexing (OFDM) is a 
method that allows transmitting high data rates over extremely degraded 
channels at a comparable low complexity. In the classical terrestrial 
broadcasting scenario, in contrast to, for example, satellite 
5 communications where we have one single direct path from transmitter to 
receiver, we have to deal with a multipath-channel as the transmitted 
signal arrives at the receiver along various paths of different length. Since 
multiple versions of the signal interfere with each other {inter symbol 
interference (IS!)) it becomes very difficult to extract the original 

10 information. The common representation of the multipath channel is the 

channel impulse response (cir) of the channel, which is the signal received 

^2 at the receiving station if a single pulse is transmitted from the transmitter. 

i if we assume a system transmitting discrete information in time 

Y: intervals 7, the critical measure concerning the multipath-channel is the 

delay Tm of the longest path with respect to the earliest path. A received 
^ ^ symbol can theoretically be influenced by Tm/T previous symbols. This 
i;:: influence has to be estimated and compensated for in the receiver, a task 
fij that may become very challenging. 

£3 Multi-path transmission of the data upon a radio channel or 

20 other communication channel introduces distortion upon the data as the 
data is actually communicated to the receiving station by a multiple 
number of paths. The data detected at the receiving station, therefore, is 
the combination of signal values of data communicated upon a plurality of 
communication paths. Intersymbol interference and Rayleigh fading 
25 causes distortion of the data. Such distortion, if not compensated for, 
prevents the accurate recovery of the transmitted data. 

Various methods are used to compensate for the distortion 
introduced in the data during its transmission upon a communication path. 

The ability to obtain reliable channel estimates affects the 
30 system performance considerably. A common way of estimating the 

channel in TDMA (time division multiple access) is to transmit a training 
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sequence and evaluate a Least square (LS) estimate of the channel at the 
receiver based on the knowledge of the training sequence. The LS 
channel estimate is basically a noisy version of the exact channel 
estimate. Hence, this technique relies on a low noise environment. 
5 Simulations show that for a uncoded system, a gap of about three dB at 
BER floor of 0.01 exists when using the LS channel estimate in 
comparison to using the exact channel estimate. This points to the 
advantages of using interpolation coefficients (with the least possible 
complexity) to enhance the LS channel estimate. 

10 The correlation properties of the channel have been used to 

enhance the LS estimate. For example in the paper authored by J. J. 
%i Vande Seek, O. Edfors, M. Sandell, S. K. Wilson, and P. O. Borjeson, "On 
J; Channel Estimation in OFDM systems," in proc. 45^*" IEEE on Vehicular 
J J Technology Conference, IL, July 1995, pp. 815-819, time correlation is 
ta used for channel estimate enhancement. A time interpolator relies on the 
^ correlation between different channel taps in the time domain, which 

requires the knowledge of the channel statistics versus time. The 
fy technique requires calculating the interpolator for every transmission 
IZ burst. The interpolator requires a matrix inversion of dimension N (the 
M size of the training sequence) for every burst which increases the system 
complexity. 

In the paper authored by J. J. Vande Beek, O. Edfors, M. 
Sandell, S. K. Wilson, and P. O. Borjeson, "OFDM Channel Estimation 
with Singular Value Decomposition," in proc. 46^*^ IEEE on Vehicular 

25 Technology Conference, Atlanta, GA, Apr. 1996, pp. 923 927, interpolation 
in the frequency domain is used to enhance the LT estimate. This 
technique suffers from increased complexity due to the requirement of a 
matrix inversion. This technique was modified to include low rank 
approximation in the interpolator to decrease complexity, however, the 

30 modified technique requires estimation of a group of dominant eigenvalues 
and eigenvectors for every transmission burst. Since performing such 
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eigendecomposition is a complex task, the modified teclinique suffers from 
complexity as well. 

In the paper authored by Y. Li, L. J. Cimini, Jr. and N. R. 
Sollenberger, "Robust Channel Estimation for OFDM Systems with Rapid 

5 Dispersive Fading Channels," IEEE Trans. On Communications, vol. 46, 
No. 7, July 1998, both the time and frequency channel statistics are used 
for interpolation. While reliance on both statistics enhances the channel 
estimate, it requires the knowledge of both time and frequency statistics 
for every transmission burst. In addition, calculations must be performed 

10 by the interpolator for every burst. Determining the channel statistics, 
every burst is also a very difficult task. This technique also requires 

^3 additional processing capacity at the receiver to estimate the channel 

J statistics from the received signal. This in turn increases the complexity of 
J the receiver. 

tS In the paper authored by Y. Li, N. Seshadri and S. 

Ariyavisitakul, "Channel Estimation for OFDM Systems with Transmitter 
ril Diversity in Mobile Wireless Channels/' /EEE JS>AC, vol. 17, No. 3, March 
I! 1999, a channel estimate for space time coding (STC) was introduced that 

5; basically evaluates the LS estimate of the channel in the time domain 
20 without doing any interpolation to avoid relying on the channel statistics. 

While the LS estimate alone without interpolation suffers from noise, in the 

presence of more than one transmitting antenna, it will also suffer from 

interference. 

In the paper authored by S. K. Wilson, R. E. Khayata and J. M. 

25 Cioffi, "16 QAM Modulation with Orthogonal Frequency Division 

Multiplexing in a Rayleigh-Fading Environment," in proc, VTC-1994, pp. 
1660-1664, Stockholm, Sweden, June 1994, a different approach for fast 
fading channels was introduced. This approach relies on adaptive 
interpolation. Use of this adaptive algorithm incurs problems related to 

30 algorithm convergence, i.e., the eigenvalue spread of the received data. 
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Such impairments as described above hinder the 
implementation of the LS channel estimator in real time applications. 

SUMMARY 

The invention presents a method an apparatus for estimating 
channels in orthogonal frequency division multiplexed (OFDM) 
communication systems. The method and apparatus allows a channel 
estimate to be determined independent of having knowledge on channel 
statistics. The method and apparatus may be implemented In OFDM 
systems having single or multiple transmitting antennas. 

In an embodiment of the invention, the method and apparatus 
is implemented in an OFDM system utilizing at least two antennas. 
Channel estimation is performed by determining and then utilizing a least 
square (LS) estimate and an interpolation coefficient for each transmitting 
antenna. According to the embodiment of the invention, the interpolation 
coefficient is determined independently from the statistics of the channel, 
i.e., v/ithout needing the channel multipath power profile (CMPP). The 
interpolation coefficient is determined by estimating the maximum delay 
encounted by the channel, calculating a maximum number of multipaths L 
by dividing the maximum delay by the transmitted symbol duration, 
creating a channel multipath power profile for the receiver using L, and 
performing a fast fourier transform (FFT) on the multipath power profile to 
generate a frequency correction vector which is used to determine an 
interpolator coefficient in the form of an interpolator matrix M. The 
interpolator matrix M is then multiplyed by an LS estimate for each 
transmitting antenna to determine the channel estimate for each channel. 

The method and apparatus provides a channel estimate, which 
is very close to the exact channel. Moreover, it can be readily applied to 
different communication systems such as MIMO (Multi Input Multi Output), 
SIMO (Single-Input Multi-Output), MISO (Multi-Input Single-Output) and 
(Single-Input Single-Output). The method and apparatus does not rely on 
knowledge of the channel statistics (either in time or frequency) to 
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enhance the LS estimate, and does not require such information. The 
interpolator is implemented mathematically by multiplying the LS estimate 
by the matrix M. 

The matrix M is required to be estimated once, hence, the 
5 technique does not require estimating M every burst and does not include 
any mathematical operation except multiplication. Consequently, the 
approach has a very limited complexity , and therefore , can be easily 
implemented. 

BRIEF DESCRIPTION OF THE DRAWINGS. 

1Q Figure 1 illustrates portions of a receiver according to an 

embodiment of the invention; 

Figure 2 illustrates portions of a channel estimator according to 
[y an embodiment of the invention; 

Figure 3 illustrates process steps performed when applying 
li interpolation according to an embodiment of the invention; 

fU Figure 4 is a flow chart illustrating process steps performed 

Q when calculating interpolation coefficients according to an embodiment of 
' the invention; and 

Figure 5 is a flow chart illustrating process steps performed 
20 when applying interpolation to estimate a channel according to an 
embodiment of the invention. 

DETAILED DESCRIPTION 

In the following description, particular embodiments of the 
invention are shown and described. A person skilled in the art will 
25 recognize that certain modifications may be made therein without 
departing from the scope and spirit of the invention as set forth and 
claimed. 
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Referring now to Fig. 1, therein is a functional block diagram 
illustrating portions of an orthogonal frequency division multiplexing 
(OFDM) receiver 100 according to an embodiment of the invention. 
Receiver ICQ includes time synchronizer 30, frequency offset corrector 32, 
5 fast fourier transform (FFT) operator 34, channel estimator 36, channel 
corrector 42, demodulator 44, deinterleaver 46, depuncturer 48, Viterbi 
decoder 50, and phase corrector 52. Phase corrector 52 includes pilot 
remover 38 and phase tracker 40. 

According to Fig. 1, a signal r(t), received over a radio channel, 

10 is input to time synchronizer 30. Time synchronizer 30 synchronizes the 
signal to the beginning of a transmission burst or block. Frequency offset 

%j corrector 32 then corrects the signal for any offset errors that occur 

J? between the transmitter local oscillator and the local oscillator of receiver 

pi! 100. The corrected signal is then input to FFT operator 34 and converted 

11 from the time domain to the frequency domain. The frequency domain 

^ signal is then input to phase corrector 52, which comprises pilot remover 

J;-; 35 and phase tracker 40. Phase correctors 52 provide an estimate of the 
ry phase to channel corrector 42. Channel estimator 36 also receives the 
IZ frequency domain signal and provides an estimate of the gain that the 
20 channel has incurred to channel corrector 42, which provides the 
corrected signal to demodulator 44. 

Demodulator 44, deinterleaver 46, depuncturer 48, and Viterbi 
decoder 50, together form the decoder function in receiver 100. 

Referring now to Fig. 2, therein are illustrated portions of 
25 channel estimator 36 of Fig, 1. Buffer 54 receives the frequency domain 
signal from FFT operator 34 and stores a training sequence from the 
frequency domain signal. A least squares (LS) channel estimate is then 
determined by performing division on the training sequence in LS 
estimator 56. Channel estimate decoupler 58 then decouples the LS 
30 channel estimate for each channel received over a separate antenna if 
more than one trasmitting antenna is being used, i.e., over each of a 
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plurality of antennas. Coefficient interpolator and channel estimator 60 
then receives each decoupled LS channel estimate from decoupler 58. 
Coefficient interpolator and channel estimator then multiplies interpolation 
coefficient for each channel by the LS estimator to obtain final channel 
5 estimates. 

To describe the functions of channel estimator 36 in the 
embodiment of Fig. 1, the case of two transmitting antennas may be used 
as an example. The embodiment however, may be implemented for any 
number N of transmitting antennas. 

10 An OFDM transmitter having two transmitting antennas (Tx1, Tx2) 

C:i transmitting to receiver 100, with receiver 100 having one receiving antenna (Rx), 

Va for a down link transmission (the general case of M transmitting antennas is 

12 straightforward) will be used in this example. Each transmitting antenna Txl , Tx2 

[y of the transmitter may use a long training sequence of length N. The training 
sequences of Tx1 and Tx2 may be represented by [ab] and [c,d] respectively, 

^ and chosen to be related as follows: 

rli ^ = ^ 

I- C = Ae'''^^ [1] 

Any number and choice of training sequences may be used. 
This description is generalized to any number and choice of the training 
20 sequences. 

The received signals for the two training sequences input to LS 
estimator 56 can be expressed as, 

=QA+7QA+ni, [2] 

z, =Q^hi-jQ^h,+n2, [3] 

25 Where is assumed to be the diagonal NxN matrix whose 

entries are the elements of A, is assumed to be the Nxl channel 
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response for the i**" (i6{1 ,2}) transmitting antenna, is assumed to be the 
Nx1 noise vector associated with the i*^ (ie{1 ,2}) received training 
sequence, and has a variance cr^. 

The least squares (LS) estimate for Tx1 and 1x2, respectively, 
output from channel estimator 58 hjand h, would be given by: 

h,, = 0.5Q, (z, + z, ) = h, + = hx + vx [4] 

h^,^ =0.5Q^(7Z2-7Z,) = h2 + -^^^^ = h2 + V2 [5] 

Where viandvz would be the new noise vectors with variance — . From 

[4] and [5], the LS estimate may be obtained by dividing the received 
training sequences with the actual ones. It can be also noted from [4] and 
[5] that the LS channel estimate is a noisy version of the exact one (i.e. 
the LS channel estimate is the exact channel response plus noise). 

According to the embodiment, the channel is estimated by 
coefficient interpolator and channel estimator 60 using a IVIMSE based 
filter to enhance the LS channel estimates represented by [4] and [5]. 
This mitigates the effect of the noise vectors in equation [4] and [5] by 
decreasing the noise energy (variance). This Is done by combining the LS 
channel estimates received from channel estimate decoupler 58 with 
suitable interpolating coefficients that are determined in coefficient 
interpolator and channel estimator 60. Mathematically, this is manifested 
by multiplying the LS channel estimate represented by equations [4] and 
[5] with an interpolating matrix M, 

h, = M.h„, f = l,2 [6] 

The MMSE interpolator coefficient M is based on the well- 
known MMSE criteria. 
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R 



= £[xy''] and x'' would be the conjugate transpose of x. 



In particular, the filter M minimizes the average error between 
the interpolated LS channel estimate h, and the exact channel response 
h, . This has the effect of preserving the useful term in equations [4] and 
[5] (i.e. h,) while minimizing the noise term (i.e. v,)- Ideally, the MMSE 
filter Mmay be written as 



/ 2 



M = R •(r+Rv,J"' = R» R+-^I t7] 



V 2 , 



Where in equation [7], it is assumed that channel responses 
corresponding to antennas Txl and Tx2 have the same correlation 
function R or equivalently the same Channel Multipath Power Profile 

(CMPP). 

The rank of R is almost equal to the number of non-zero taps 
in the CMPP, which is usually less than the overall dimension N, and-the 
entries of R represent the correlation between the different components 
of he,J =U, the more correlation between carriers we have, the more 
enhancements we expect from the interpolator. In a typical OFDM system 
there is a correlation coefficient of about 0.9 between each two adjacent 
carriers. 

The following algorithm can be used to interpolate the channel 
if the channel statistics manifested in CMPP is known: 

\npuV. h,;,,z=l,2. 
Output: h,w = l,2. 
Algorithm: 

For a particular radio channel knowing CMPP, find 

R = Toeplitz[FFT{CMPP)] . 
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Knowing the noise variance, substitute in [7] to get M. 
Substitute in equation [6] to get ji,.,? = l,2. 

It is to be noted tliat the CMPP is not available at the receiver. 
Hence, the above algorithm is replaced by an algorithm according to the 
method and apparatus of the invention. 

It appears clear from the analysis of [7] that the interpolator 
depends on the channel correlation function R . R is the Toeplitz matrix 

built from the FFT of the CMPP, consequently the solution will depend on 
the channel multipath power profile (i.e. CMPP). 

The embodiment of the invention provides an approach that 
almost does the same job as the exact MMSE interpolator without 
depending on the knowledge of CMPP (or equivalent the channel 
statistics) at the receiver. According to the embodiment, the above 
algorithm is replaced by an algorithm that may be performed independent 
of knowledge of the CMPP. The following Lemma may be used to 
describe the method and apparatus. 

Lemma: 

If H, = IDFTiii), i = 1,2 , H, ^ = IDFT(h, ,^ ), i = 1,2, a is the vector 
constructing the teoplitz matrix R (the first column in R) and 
(p^ik) = {lDFT(a)),,k = 1,2,---,N then equation [6] corresponds in the time 

domain to 

a = 'i'*H,, [8] 

0 

0 TO) ••• 

Where . . 

0 0 •• 



0 
0 

T(A^) 



and 
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Proof: 



The expression in [8] can be proved by recalling from [4] and 

[5] that, 

h,,, =h, + v,,J=l,2 [11] 

Applying the IDFT operator to [11] we get, 
H,,, =:a + V,,i=l,2 [12] 

where h, = IDFT(hX i = 1,2 and due to the orthogonality of the 
/DFT operator, the new noise components are also independently 

identically distributed (iid) but with a covariance matrix —I. Solving for 

the MIVISE filter F that estimates H,from H,,Jn equation [12], we get, 

F = R„„ •R"' [13] 



where R 



R 









0 


0 




0 


0 



0 
0 



[14] 



The expression of R^^ results from the fact that the channel 

15 coefficients are uncorrected for different paths, hence the off-diagonal 
entries in R^,^^ vanish or equivalently, R^,^^ is a diagonal matrix. The 

diagonal entries represent the power in each path, i.e. the components of 
the CMPP. Substituting equation [14] in equation [13], then equation [8] 
follows, 

20 Equation [8] indicates that the function of the interpolator is 

equivalent in the time domain to scaling the component of the LS 
channel estimate for each transmitting antenna with ^{k). The person 
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skilled in the art will recognize that the number of multipaths in the 
channel is usually much less than the number of carriers N. Hence, only 
few taps of the LS channel estimate in the time domain are carrying useful 
energy while, the rest are only noise. Stated differently, referring to 
equation [12], the useful term in equation [12], u,, has few nonzero 

entries while the entries of the noise term Vi are all nonzero. Since W{k) 
and a have nonzero entries at the same positions, scaling the 

component of the LS channel estimate with T(^) basically preserves 
the useful part in equation [12] (i.e. hJ and eliminates a major portion of 
the noise part (i.e. vj- Based on this, it can be noted that: 

Since the value of the non-zero '¥(k) in equation [8] is close to 

one (even at very low SNR value as — ««=';.(^)), then the exact value of 

the rnuitipath profile used at the receiver is irrelevant and what really 
matters is the positions of these taps. In other words, we can achieve 
almost the same performance if the receiver used a Receiver Multipath 
Power Profile (RMPP) that differs from the channel one (CMPP) as long as 
it does not miss a tap in CMPP (i.e. as long as there is no zero entry in 
RMPP which corresponds to a nonzero entry in CMPP). 

f the receiver misses a tap that exists in the channel than it is 
scaling some received path by a zero value or equivalently eliminating 
some of the received energy. It is to be expected that such a scenario 
would deteriorate the interpolator performance. 

If the receiver does not miss a tap in the channel, however, it 
adds more taps than those really exists, it is basically collecting noise at 
these taps. Simulations show that the influence of picking up such noise is 
not significant since Lch « N 

The maximum number of channel taps L^^^hat can exist is so 

well defined, i.e. the ratio between the channel multipath spread Tm and 
the symbol duration T. Thus, a scenario that achieves most of the 
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interpolator performance with mucin less complexity is to fix a multipath 
power profile at the receiver that basically includes a number of taps equal 
to L^, . In such case, the RMPP will never miss a tap that is in CMPP. 

Based on the knowledge of L^^ , the coefficient interpolator and 
5 channel estimator 60 will use a RMPP covering all the expected taps In 
CMPP. The values of the interpolation coefficients can then be 
determined (based on only knowing 4J. The coefficient interpolator and 
channel estimator 60 then would use these coefficients to interpolate the 
LS channel estimate. It is to be noted again that the same coefficients are 
10 to be used every burst, so the coefficient interpolator and channel 

estimator 60 need not to calculate M (and hence find the inverse of NxN 
'if matrix) every burst. 

fl According to the embodiment, when a RMPP that consists of 

4a taps is chosen with any power values. R - FFT{RMPP) Is then used in 
l5 the algorithm instead of R. 

5=}1 Referring now to Fig. 3, therein are illustrated the process 

Si steps when calculating interpolation coefficients according to an 

ci embodiment of the invention. A received time signal consisting of the 

training signal is convoluted with the channel plus White Gaussian Noise 
20 (WGN) (1). The time signal is then converted to the frequency domain via 
FFT operation (2) in FFT operator 34. The LS estimator 56 multiplies the 
received signal in the frequency domain by the conjugate of the training 
sequence (3) to result In a noisy version of the channel response. 
Coefficient interpolator and channel estimator 60 takes the LS estimate in 
25 the time domain (4). Due to scaling performed according to equation [8], 
the coefficient interpolator and channel estimator 60 scales the first 
4^ components using ones and It replaces the last N- components by 
zeros (5). This process has the effect of suppressing a lot of noise 
components while not affecting all the channel components since the 
30 channel can only exist at some positions in the first L^^ components. The 
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new (less-noisy) estimate is then transformed to the frequency domain (6). 
Consequently, the interpolator acts as a low-pass filter but in the time 
domain. 

Referring now to Fig. 4, therein is a flow chart illustrating 
5 process steps when calculating the interpolation coefficient according to 
an embodiment of the invention. As already mentioned, it will not be 
necessary that a calculation be performed every burst but instead it can 
be done once as long as the channel multipath spread Tm is constant. 
The multipath spread Tm for those channels is pre-known to the designer 

10 usually from intensive measurements that had been done on such 
channels. Hence, the requirement of knowing Tm adds no burden to the 
receiver complexity. 

In block (10) an estimate of the maximum delay encountered 
|l| by the channel is performed. From block (10) the maximum number of 

11 multipaths L can be calculated by dividing the maximum delay 
encountered by the channel Tm by the symbol duration T (12). In block 
(14), a receiver multipath power profile is created. Next, in block (16) by 

!^ performing an FFT operation on the receiver multipath power profile, the 
ij frequency correlation vector is found. Next, in block (18), the interpolator 
matrix M is calculated by constructing the teoplitz of \|/. 

If M is multiplied by the least square channel matrix obtained 
by the process described in Fig. 5 the final estimate of the channel is 
obtained. 

Referring now to Fig. 5, therein is a flow chart illustrating 
25 process steps when applying interpolation according to an embodiment of 
the invention. The process described in figure 6 is a burst by burst 
process to obtain the least square channel estimate. The received signal 
r(t) is put into the frequency domain by the FFT operation (20) and the 
training sequence is extracted from the preamble of the burst (22). A least 
30 square channel estimate is obtained by dividing the received training 

sequence by the exact training sequence (24). Block (26) exists only in the 
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case of multiple antennas case and comprises the step of decoupling the 
different channels corresponding to the different transmitting antennas. 

In block (28) a complex matrix -vector multiplication is 
performed, by multiplying the least square channel estimates and the 
interpolating coefficients to estimate each channel. 

Thereby, a manner is provided by which to communicate data 
on a channel susceptible to distortion. When utilized, an improved and 
simplified communication method of communications is permitted. The 
preferred descriptions are of preferred examples for implementing the 
invention, and the scope of the invention should not necessarily be limited 
by this description. 
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